SSC98-1X-6 


Integrated Power and Attitude Control System (IPACS) 

Ted Michaelis Code 573 - , - j/ i 

NASA-GSFC 
Greenbelt, MD 2077 1 
301-286-5939 

theodore.d.michaelis. l@gsfc.nasa.gov 

Abstract Recent advances in materials, circuit integration, and power switching have given the concept of dynamic 
energy and momentum storage important weight, size, and operational advantages over the conventional momentum 
wheel-battery configuration. Simultaneous momentum and energy storage for a three axes stabilized spacecraft can 
be accomplished with a topology of at least four wheels where energy (a scalar) is stored or retrieved in such a 
manner as to keep the momentum vector invariant. This study, instead, considers the case of two counter-rotating 
wheels in one axis to more effectively portray the principles involved. General scalable system design equations are 
derived which demonstrate the role of momentum storage when combined with energy storage. 


Introduction 

General scalable system design equations are derived 
that demonstrate the role of momentum storage when 
combined with energy storage. Scaling to “Smaller 
than Small”- the motif of this conference - may 
require innovative design and a search for suitable 
vendors. 

Rotor and Enclosure Mass 

Parameters are defined in the glossary and units are 
SI. Two areas of mass minimization are considered. 
One is the derivation of the expression for minimum 
mass sum of rotor and enclosure masses. It is 
assumed a pressure enclosure is required 

corresponding to an unmanned application. 
Additional mass for catastrophic containment is not 
considered, A second area is mass minimization of 
the motor-generator (M/G) and that solar array mass 
fraction required to support its copper loss. 

Mass comparisons are made between a conventional 
momentum wheel-battery system, battery 

replacement by dual wheels, and an IPACS system. 

Practical implementation is discussed and EPACS 
advantages listed. 

Configuration 

A pair of counter-rotating wheels on three orthogonal 
axes is simplest to perceive. Another arrangement 
utilizes large angle control moment gyros with 
important power and weight advantages. A third 
geometry employs four wheels with momentum 
vectors perpendicular to the faces of a tetrahedron. 
This study is largely confined to the study of two 
counter-rotating wheels. Extension to other 
geometries is straightforward 


The rotor is assumed to be a hollow cylinder with 
outer radius r , . The height and inner radius are 

defined in Table 1. With k x and k 2 fixed, rotor 
shape remains invariant as size changes. Rotor inertia 
and mass are given in Eqs. (1) and (2). 



(1) 

M = npk 1 (1 - k] )r, 3 

(2) 


The moment of inertia required by Eq. (4) is 
obtained with the aid of Fig. 1. and the fundamental 

Eq. (3). 



This gives the outer radius tangential velocity and for 
the representative values of tu and p given in the 

table yields 1643 meters/sec. or 1.643 Km/sec. 
(3700 miles/hr.). <j is the operating stress which 
may be smaller for a given application than the value 
given in Table 2. It should be as high as possible 
consistent with safety and operating life. 



Two Counter-rotating Wheels 



r,E + {\+k)H 


J = -±- 


(1 ~k 2 )~ 


(4) 


Eq. (4) is the inertia required to support H and E with 
angular velocity ranging between kco s and (t ) s . 

This is set equal to Eq. (1) after substitution of Eq. 

(2) into Eq. (1). A power equation inr t results: 

npk^X-k^l-k 1 )— r* -2Er x -2(1 + *)// I— = 0 
P \P 


For the values given in Table2, r, = 0.07749 and 
the rotor weight has decreased from 2.52 Kg to 1.91 1 
Kg. But k 2 has increased from 2 to 3.97. 


Motor-Generator Size 

The motor-generator (MG) size is based upon its 
copper loss. The I 2 R dissipation is invariant with 
winding impedance for a given size but decreases as 
motor size (and mass) increase. This power loss must 
be provided by the solar array with its attendant 
mass. The goal is to find an analytical expression that 
minimizes the sum of these two masses. 


Eq. (5) 

Given the two rotor shape parameters ,k x and k 2 , 
Eq. (4) uniquely determines the outer radius. There is 
no mass minimization. With the values in Table 1, 

r x =0.0842 m. and the rotor mass, from Eq. (2), is 
2.52 Kg. 


The basic steady state DC motor equation is: 

V=IR +k B m (8) 

Multiplying by I; 

VI=P =I 2 R + k B Ico (9) 


Another approach is to let k 2 be a free parameter 
and consider a housing weight constraint. Then a 
radius can be found that minimizes the sum of rotor 

and housing masses and also fixes k 2 . Its value, 
however, cannot be too large, leading to a rod shape, 
difficult to dynamically stabilize. The governing 
equation is: 


2[ E + (\ + k)y\- ] 

M + M„= r ±J_P_ + (6) 

(i-* 2 Xi + *, 2 )- 

P 

The first term is derived fromEqs. (1), (2), and (4) 
and the second term results from analysis of a 
working momentum wheel. Setting the derivative of 

Eq. (6) with respect to r x equal to zero, the value 
which produces minimum mass is: 


2 H 


3* 3 (l-*)(l + *, 2 ) 1 


(7) 


Substitute for I using T =k T l where k T is the torque 
constant. In SI Units k B = k T = k i 


T 1 

P= — -R + Tco (10) 

*< 2 

This is generalized by invoking the motor constant 
which is commonly available from the vendor. 



( 11 ) 


I >=-^ + 7© (12) 

A* 

The first term is the copper loss, independent of 
winding impedance, and the second term is the 
mechanical power. It remains to find the relationship 

between K M and motor mass, M m . 


Let b= — — (13) 

A survey of 29 Inland motors ranging in mass from 
45 grams to 8.3 Kg gave a value for b of 5.23 with 
a standard deviation of 1.28 with rare earth magnet 



motors having smaller values in this set. Finally the 
motor copper loss is: 


P = 


{ Tb V 


M r 


Eq. (14) is solved for M m : 


(14) 


M m 


Tb 

V? 


(15) 


Comparison s 


Three cases are examined. 

a) Energy storage only 

b) Conventional battery-wheel system 

c) IPACS 

Energy Storage On ly 

In Eq. (4) H=0 and then fromEqs. ( 1 ) and (4): 


The solar array has a figure of merit; c watts/Kg.. 
Hence: 


p 2 E 

M= — T ~Z 

a (1 + ^X1-* ) 


(19) 


M p 



(16) 


This is the incremental array mass required to support 
the MG copper loss. The sum of Eqs. (15) and (16) 
is differentiated with respect to P, set to zero, and the 
resulting value of P given by Eq. (17) is substituted 
into Eq. (15). 


P = 


TbcX i 

< 2 J 


(17) 




l 

2b 2 T 2 X 

c 


(18) 


This is independent of r, and k 2 . In the limit as 
k -> 0 and k x —> 1 

M =£-E (20) 

<J 

Eq. (20) represents the ideal minimum mass thin 
shell rotor., useful for comparison. 

Using the values in Table 2 Eq.. (1 provides a mass 
of 

0.682 Kg. for the mass of each of two counter- 
rotating wheels. 

Table 1 

Wheel Parameters for 300 WHrs 


b and c are fixed so that the MG mass increases as 
the 2/3 power of torque . To put this in perspective, 
the benchmark often used here, is 300WHrs storage. 
The discharge power is then 900 watts assuming30 
min. eclipse and 90 minute orbit. If rotor speed is 
4000 radians/sec. (half full speed), then the torque is 
0.225 Nm and with b given above and c=12 watts/Kg 
the MG mass given by Eq. (18) is 0.61 Kg. From Eq. 
(17) the power is 3.68 watts. 

It might seem that increasing the wheel depth of 
discharge would decrease rotor mass. However this 
would increase motor mass since, for a constant 
power withdrawal during eclipse, (P=Tft> ) the lower 
the speed the greater the torque. Another negative 
effect is the increased complexity of the power 
conditioning circuits required to utilize power over a 
greater voltage range as well as the array element 
complexity needed to efficiently charge at a lower 
voltage. 


Item 

Mass Kg 

rotor 

0..68 

MG 

0.61 

bearing 

0.55 

housing 

0.62 

Total 

2.46 


For two wheels the mass is 4.92 Kg. Energy storage 
wheels are not temperature sensitive as batteries are, 
have no limit on depth of discharge as a function of 
number of cycles and no time dependent failure 
modes.. It would be advantageous on GEO spacecraft 
since it could be completely de-energized except at 
equinoxes. 




A state of the art nickel-hydrogen battery is used for 
comparison. It will be de-rated to 10 year LEO use 
where on each discharge cycle it is discharged 
to only 67% of its rated capacity. A nickel-hydrogen 
battery rated at 15 AH at 28 volts weighs 8.9 Kg. A 
30 ah battery, scaled up to 32.1 AH and de-rated by 
2/3 will provide 300 WHrs usable energy and weigh 

19.0 Kg. The dynamic energy storage mass is 
therefore about 25% of battery mass. Array element 
switching and battery unique circuit masses are not 
considered; both are small and have similar values. 
The mass saving for direct battery replacement is 

14.1 Kg. 


Conventional Battery- Wheel System 

A mass estimate for a conventional battery and 
wheel system is made to provide a comparison with 
an IP ACS system. The battery mass will be the same 
as calculated in the previous section (19.0 Kg). A 
three axis stabilized system will require a 50Nms 
wheel on each axis with the following mass. 1,2 


M ^ “ 2.8// 8 = 12.1 Kg (21) 


Spin Axis Orientation Schem e 

Magnetic bearings have the design latitude to provide 
a small degree of freedom, typically+_ 1 degree in 
spin axis orientation. An interesting but limited 
system may be designed on this basis. A 
tetrahedronal array of magnetically suspended 
wheels is considered where each has a small 
transverse angular motion capability, a . All four 
wheels operate at the same speed over the range from 

ko) s t oO) s With nominal orientation the net 

momentum is zero. To keep it simple and 
symmetrical, momentum is desired in the direction 
of one of the spin axes. This net momentum is 
obtained by symmetrically changing the spin axis 
orientation of the three remaining wheels. This is 
probably a worst case as well. Wheel speeds are 
assumed at the discharged value, which is worst case. 
Only the finally derived equation is presented; based 
upon some of the foregoing relationships. 


H = -kJip\ 


E* 

5 

nk 2 (\-k*)(\-k 2 )(7~ 2 


4i 

a(l-^a) (22) 


The total mass for this system is then 55.3 Kg. 
(Other ACS masses common to all systems, such as 
momentum de-saturation torquer bars are not 
included.) 

IPACS 


Two wheels for each of three axes is base-lined for 
direct comparison with the conventional system. It is 
assumed that 50 Nms capability is required on each 
axis and that 300 WHr storage is divided between 
three axes for 100 WHr per axis. Note that E in the 
equations has the units of joules. Other parameters 
are listed in Table 2. Then from Eq. (6) the rotor 
radius is 0.07749 and with this value used in Eq. (7) 
the motor and housing mass equal 1.228. From the 
argument following Eq.(18), the torque is one third 
the former value. Then Eq. (16) givesan MG mass of 
0.295. Including the bearing estimate from Table 1 
the single wheel mass is 2.07 Kg, and for six wheels 
12.4 Kg. This is 23% of the conventional system 
mass. 

One reason for the advantage is poor battery 
utilization (only a third of its capacity can be used). 
Another is the efficient use of the high performance 
inertial rotors. The mass saving over a conventional 
system is 43 Kg. 


H depends upon energy, carbon fiber properties, rotor 
form factors, and per unit speed at discharge. The 
curious 4/3 power relationship between H and E is 
due to the increase in moment of inertia and outer 
radius with E.. With values from Table 2, N is 3.6 
Nms for E= 5E6 joules. It seems that bigger is better 
in this instance but a parametric analysis of 
disturbance momentum as energy (and presumably 
spacecraft size) increase, is needed. 

Rotor Shap e 

There are differing opinions about rotor shape with 
advocates of rod shapes, disk shapes, and a compact 
shape with height equal to outer diameter. The 
foregoing equations show how these parameters 
influence the design on a theoretical basis. More 
important is choosing viable mechanical designs that 
produce a rotor that stays together with minimal 
modal problems. Then the system response to 

successfiil rotor geometries can be assessed via these 
equations. Design need not be confined to a 
cylindrical section. For those containing a central 
shaft the cross sectional shape approaches a truncated 
Gaussian curve for most efficient use of fiber. 



Conclusions 


Simple battery replacement by two counter-rotating 
energy wheels saves appreciable mass but the IP ACS 
system provides the greatest mass saving, even with 
large momentum storage capability. Coupled with 
greater transient discharge capability and temperature 
insensitivity, the energy wheel will be a major future 
contender. 

Another advantage of IPACS is the elimination of 
dynamic unbalance torque by virtue of the magnetic 
bearings which can be designed to have low gain at 
higher frequencies such that the rotor will rotate 
about its mass center rather than geometric center. 

The intent of this paper is to provide a primer on 
design of an energy wheel system in general scalable 
terms suitable for large and small systems. As such, 
the numerical value of carbon filament modulus, for 
example, is not to be taken as authoritative. 
Numerical values are illustrative and in the ball park. 
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Table2 

Symbols 


Symbol 

Description 

Units 

Value 



ATgV wart 


b 

MG figure of merit 

Nm 

5.23 

c 

Solar array figure of merit 

watts 

Kg 

12 

E 

Stored energy capacity 

joules 



Example 

joules 

1.08E6 

h 

Rotor height 

m 

- 

H 

momentum 

Nms 

50 

J 

Single rotor moment of inertia 

Kg- m 1 


k 

Per unit lower speed limit 

- 

- 


Example 

- 

0.5 

ki 

Ratio if inner rotor radius to outer radius 

- 

- 


Example 

- 

0.75 

k 2 

Ratio of rotor height to outer radius 

- 

- 


Example 

- 

2.0 

k 3 

Housing mass penalty 

Kg 



Example 

m' 

480 

k4 

MG back emf and torque constant 

Kg- m 1 
As 3 

' 

kr 

MG torque constant 

Nm! Amp 

- 

k B 

MG back emf constant 

- 

Kg 

Km 

Kilograms 
Motor constant 

volt - sec 

Kg 

Nmf^watt 

; 

M 

Rotor mass 

Kg 

- 

M h 

Housing mass 

Kg 

- 

M m 

MG mass 

Kg 

- 

r i 

Rotor outer radius 

m 

- 

r 2 

Rotor inner radius 

m 

- 

CT 

Carbon fiber maximum working modulus 

Pascal 

4.139E9 

P 

Carbon fiber mass density 

Kg/m* 

1533 

(0 

Rotor radian velocity 

1/s 

- 

CO, 

Rotor maximum radian velocity 

1/s 

- 

a 

Tetrahedron array; spin axis orientation control 

radians 

K 

180 




